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i-C Abstract 

O 



Classical Heisenberg antiferromagnets with uniaxial exchange anisotropy and a cu- 
bic anisotropy term in a field on simple cubic lattices are studied with the help 
of ground state considerations and extensive Monte Carlo simulations. Especially, 
we analyze the role of non-collinear structures of biconical type occurring in addi- 
tion to the well-known antiferromagnetic and spin-flop structures. Pertinent phase 
diagrams are determined, and compared to previous findings. 
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1 Introduction 



Uniaxially anisotropic Heisenberg antiferromagnets in a magnetic field have 
been studied quite extensively in the past, both experimentally and theoreti- 
cally. Typically, they display, at low temperatures, the antiferromagnetic phase 
and, when increasing the field, the spin-flop phase [I]. A prototypical model 
describing these phases is the Heisenberg model with a uniaxial exchange 
anisotropy, the XXZ model 

Hxxz = Jj2[A(S*S* + SfS?)+S*S*} - H^St (1) 



where J is the exchange coupling between classical spins, (Sf^, S^-p ^i(j))j °f 
length one at neighboring sites, i and j, of a simple cubic lattice, A is the 
exchange anisotropy, 1 > A > 0, and H is the applied magnetic field along 
the easy axis, the z-axis. 
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The phase diagram of the model has been studied several years ago, using 
mean-field theory [2], Monte Carlo simulations [3], and high temperature se- 
ries expansions [4], suggesting the transition between the antiferromagnetic 
(AF) and spin-flop (SF) phases to be of first order and the boundaries of the 
paramagnetic phase to the AF and SF phases to be continuous transitions in 
the Ising and XY universality classes. Based on renormalization group analy- 
sis in one loop order, a bicritical point in the Heisenberg universality class has 
been proposed, at which the three different phases meet |5|6| . However, this 
scenario has been scrutinized when doing renormalization group calculations 
in high loop order [7j, where the bicritical point is found to be unstable against 
a 'tetracritical biconical point' [6], which, in turn, may be unstable towards 
transitions of first order in the vicinity of the multicritical point of the three 
phases. The seemingly conflicting descriptions may be reconciled by a new 
renormalization group analysis in two loop order [5]. 

As has been noted very recently [5)10)11] , not only AF and SF phases, but 
also 'biconical' (BC) structures, see Fig. 1, may play an important role in the 
XXZ model. Indeed, such BC structures are degenerate ground states at the 
critical field separating AF and SF configurations at zero temperature. For 
the XXZ model on a square lattice, these degenerate BC fluctuations seem to 
lead to a narrow disordered phase intervening between the AF and SF phases 
at low temperatures [9lll2|ll3] . The importance of BC structures for the three- 
dimensional XXZ antiferromagnet, where they are also present as degenerate 
ground states at the special field, had not been studied in any detail so far. 

Biconical structures may be ground states even in a finite range of magnetic 
fields, giving then rise to an ordered BC phase at low temperatures, when intro- 
ducing in the XXZ model further anisotropy terms or longer-range exchange 
interactions, as it is known for many years [T31IT5] . This feature has been con- 
firmed in recent simulations when adding a quadratic single-ion anisotropy to 
the XXZ model on the square lattice [TUfTTj . The related lowest-order single- 



(a) 





Fig. 1. Spin orientations on neighboring sites showing antiferromagnetic (a), bicon- 
ical (b), and spin-flop (c) ground state structures of the XXZ model. 



ion term of cubic symmetry may be written in the form [T6l[T 



12\4 



(2) 



where F denotes the strength of the 'cubic anisotropy'. The sign of F de- 
termines whether the spins tend to align along the cubic axes, for F < 0, 
or, for F > 0, in the diagonal directions of the lattice. Because of these ten- 
dencies, the BC structures show no full rotational invariance in the xy-plane 
perpendicular to the easy axis, the z-axis, in contrast to the XXZ case. 

When an ordered BC phase exists at low temperatures, for example due to 
the cubic anisotropy, intricate multicritical behavior may show up, including 
a tetracritical biconical point, at which the AF, SF, BC, and paramagnetic 
phases meet, as has been discussed before [Toll8lll9ll2U] . applying mean-field 
theory and renormalization group arguments. 

Experimentally, many antiferromagnets with uniaxial anisotropy have been in- 
vestigated, quasi two-dimensional magnets [21|22| 25] as well as three-dimensional 
magnets such as GdA10 3 , NiCl 2 4H 2 0, MnF 2 or Mn 2 (Si,Ge)S 4 [23I25I25I27I2S] . 
While we shall deal here with the theoretical analysis of the models, results 
may turn out to be useful for interpreting specific experiments in future work. 

The aim of our paper is to study, especially, the role of biconical, non-collinear 
structures in three-dimensional classical Heisenberg antiferromagnets with 
uniaxial exchange anisotropy and cubic anisotropy. Both ground state con- 
siderations and Monte Carlo techniques are applied. 
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Fig. 2. Ground states in the (F, H)- plane for the full Hamiltonian, 7~Cf= TCxxz + 
TCca, with exchange anisotropy A= 0.8. The full circle denotes the highly degenerate 
point in the XXZ model. 



The paper is organized as follows: First, results of ground state calculations 
are presented, and then Monte Carlo findings on thermal properties and phase 
diagrams will be discussed, for the XXZ model and its extension. The article 
will be concluded by a short summary. 



2 Ground state properties 



The ground states of the full Hamiltonian, Ti{= Hxxz + ^ca, eqs. (1) and (2), 
may be determined by minimization of the energy, e.g., with respect to the 
azimuthal angle, 0, i.e. the angle between the projection of the spin vector in 
the xy-pl&ne and the x-axis, and the z-component of the spin vector [T¥P^] . 
As usual, we consider spin structures with two sublattices, A and B, where 
neighboring sites on the cubic lattice belong to different sublattices. In general, 
the minimization may be easily done numerically. 

In Fig. 2, resulting ground states in the (F, H )-plane are depicted, setting the 
exchange anisotropy A equal to 0.8, as before PI91TT2"] . 

In the case of a vanishing cubic anisotropy, F = 0, i.e. for the XXZ antiferro- 
magnet, one encounters, when increasing the field, AF, SF, and ferromagnetic 
ground state configurations. In complete analogy to the model on a square 
lattice, the ground state is highly degenerate in BC structures at the critical 
field H c separating the AF and SF structures, see Figs. 1 and 2. For cubic 
lattices, one has H c = 6Jy/l — A 2 . The non-collinear BC spin structures may 
be characterized by two different polar or tilt angles of the spin vectors on the 
two sublattices A and B, 9^ and 9^, as shown in Fig. 1(b). Obviously, for the 
SF structures, the angles are identical, 9^ = 9# = Qsf, depicted in Fig. 1(c). 
At H c , the two tilt angles of the degenerate BC structures are interrelated by 



arccos 




(3) 



Fig. 3. Sketch of discretized biconical structures for antiferromagnets with cubic 
anisotropy, H{, in the case of (a) F > 0, BC2, and (b) F < 0, BC1. 



interpolating continuously between the AF and SF configurations, when vary- 
ing one of the tilt angles [9lfT0] . Of course, the rotational invariance of the spin 
components in the xy-pl&ne leads to an additional degeneracy in the BC and 
SF configurations. 



Let us now consider positive values of F, F > 0. For vanishing field, H = 0, 
and F/J > 0.3, see Fig. 2, the simple AF structure is replaced by a tilted anti- 
ferromagnetic (TAF) configuration. Spins on the two sublattices point still in 
opposite directions, but not along the easy axis. Indeed, the spins are oriented 
more and more towards diagonals of the lattice with increasing positive cubic 
anisotropy F |29| . Applying, for F > 0, an external field, BC structures may 
occur, in which neighboring spins have different tilt angles, 0,4 and ©#. In 
contrast to the XXZ case, there is no full rotational symmetry in the azimuthal 
angles, 0, of the spins. Instead, for F > 0, the spins order along the diagonals 
of the lattice, and is discretized, taking the values, (2n + l)7r/4, where n is 
an integer. The xy-spin components are ordered antiferromagnetically, 4>b= 
4>a + 7T. The resulting discretized biconical structure, BC2, is sketched in Fig. 
3(a). Obviously, the SF configurations are discretized in the azimuthal angle 
in the same way as the BC structures. When fixing F and varying the field 
H, the two tilt angles, ©^ and 0#, change continuosly when going from the 
AF (or TAF) to the BC2 structures. In the BC2 region, the tilt angles, Q A 
and Qb, vary continuously as well. On the other hand, they seem to jump 
when going from the BC2 to the SF region. Note that the tilt angle of the 
SF structure, Qsf, at the border to the BC2 region, see Fig. 2, depends only 
weakly on the strength F of the cubic anisotropy. 

For F < 0, the azimuthal angles of the spins, both for biconical and SF 
structures, take the values = nii/2, n being an integer, reflecting the fact 
that the cubic anisotropy now favors alignment of the spins along the cubic 
axes. The xy-spin components are ordered antiferromagnetically, <Pb= 4>a + ^, 
as for F > 0. Biconical ground states occur for sufficiently strong negative 
values of F, as shown in Fig. 2. The resulting BC1 structure is shown in Fig. 
3(b). The tilt angles ©^ and 0# change continuously in the BC1 region, when 
varying the cubic anisotropy F and the field H . At the boundaries of the BC1 
region to the SF, AF, and ferromagnetic regions the tilt angles seem to jump, 
except between the BC1 and SF regions near the triple point of these two and 
the ferromagnetic regions, see Fig. 2. 



3 Monte Carlo simulations 



In our simulations of thermal properties of the XXZ and the full model, 7^f, we 
applied the Metropolis algorithm with single-spin-flips. Lattices of L 3 sites, 
employing periodic boundary conditions, with L ranging from 4 to 32, allowed 



to do useful finite-size analyses. To obtain thermal averages, we performed, 
at fixed model parameters and temperature, several independent runs, with 
distinct random numbers, each run consisting of, at least, 10 7 Monte Carlo 
steps per site for the larger systems. 



We recorded standard thermodynamic quantities such as the specific heat, 
C, the (absolute) total magnetization, m, longitudinal, m z st , and transverse, 
staggered magnetizations, being the order parameters in the AF and 



m 



st 1 



SF phases, as well as the corresponding susceptibilities, Xi Xsv an d Xst ■ in 
addition, we computed the Binder cumulants [30] of the two different order 
parameters, U z and U xy , and related histograms. To gain microscopic insights, 
especially, on BC structures, we also recorded probability functions of the tilt 
angles, such as the probability P2(@a, ©b) for finding the two angles, 0a and 
0#, at neighboring sites, and the probability p(Q) for encountering the tilt 
angle 6 



Let us first consider the XXZ model on a cubic lattice, eq. (1). Its phase 
diagram has been determined, using Monte Carlo techniques, already some 
time ago [3]. Certainly, in present simulations the accuracy has been improved 
significantly. The resulting phase diagram is depicted in Fig. 4, where the phase 
boundaries have been estimated by standard finite-size extrapolations [3T] for 
thermodynamic quantities and the Binder cumulant. The phase boundary lines 
deviate somewhat from the ones estimated before [3], which we attribute to 
the improved statistics of the present study. 



We also identified the type of transitions and, in the case of continuous transi- 
tions, their universality classes, by monitoring, e.g., the finite size-dependences 
of the peak heights of the susceptibilities and specific heat, to estimate critical 
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Fig. 4. Phase diagram of the XXZ model on a cubic lattice with exchange anisotropy 
A = 0.8, as obtained from Monte Carlo simulations described in the text. 



exponents, see below. 



In agreement with previous suggestions [3|T3] . the transition between the AF 
and SF phases is found to be of first order, with, for instance, the maxi- 
mal staggered susceptibilities increasing with system size proportionally to 
L 3 , characteristic for such transitions [32J. Moreover, as a remainder of the 
degeneracy of the ground state at H c , biconical structures prevail near that 
transition at low temperatures, as one may easily observe in P2, see Fig. 5. 
As included in the figure, the dominant BC fluctuations are quite close to 
those expected from the degenerate ground states, eq. (3). Note that p 2 dis- 
plays simultaneously local maxima at positions belonging to the AF and SF 
structures, as shown in Fig. 5. The phenomenon gets more pronounced when 
increasing the system size, and it reflects coexistence of the AF and SF phases 
at a first-order transition. This behavior is in marked contrast to that of the 
XXZ antiferromagnet on a square lattice, where the BC fluctuations seem to 
lead to a disordered phase between the AF and SF phases 



The transition from the paramagnetic phase to the AF phase is found to be 
in the Ising universality class, while the transition from the paramagnetic 
phase to the SF phase is found to belong to the XY universality class. For 
instance, we estimated canonical [3T] asymptotic critical exponents describ- 
ing the size dependence of the peak height of the staggered susceptibilties, 
(7/V) 2 = d(ln x z ) / d(ln L) , L —>■ oo, and analogously for the transverse sus- 
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Fig. 5. Grayscale representation of the probability p2 for finding the tilt angles Qa 
and ®b at neighboring sites in the XXZ model with A = 0.8 in the vicinity of the 
transition between the AF and SF phases, at ksT/J = 0.7 and H/J = 3.6925. The 
solid line corresponds to eq. (3). The system size is L = 24. 



ceptibility, see Fig. 6. Our estimates agree nicely with the known, accurate 
values obtained from renormalization group calculations in high loop order 
for both universality classes [33] • Similarity, our estimates [29J for the crit- 
ical Binder cumulants U* and U* at the two types of transitions are close 
to those obtained for the three-dimensional nearest-neighbor Ising and XY 
models [34.35]. Note, that the critical Binder cumulant may change, within a 
given universality class, due to, e.g., spatially anisotropic interactions [36|37j . 
This possible pitfall seems to play no role in the XXZ model. 



Analyzing critical exponents and critical Binder cumulants, we locate the mul- 
ticritical point at k B T/J = 1.025 ± 0.015 and H/J = 3.90 ± 0.03, improving 
the previous estimate [3]. To obtain the estimate, we proceed as exemplified 
for xlt i n Fig. 6. We determine the effective, size-dependent critical expo- 
nent ( i y/u)lff(L) = d(\n x z ) / d(\n L) of the height of the peak in the staggered 
longitudinal susceptibility xlt f° r the simulated system sizes at various fixed 
temperatures close to the boundary of the AF phase. One observes a pro- 
nounced increase of the effective exponent in a small interval of temperatures. 
This behavior signals the change from a continuous transition of Ising type, 
with (j/v) z ~ 1.9635 [33], to a transition of first order, where the asymptotic 
exponent is expected to be 3. A pronounced increase occurs simultaneously, 
within the temperature resolution of the simulations, for the effective expo- 
nent of the staggered transverse susceptibility [29]. Accordingly, it indicates 
a transition of first order between the AF and SF phases. The ratio of the 
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Fig. 6. Effective critical exponent {^/v) z e rAL) = d\n\ z /d\nL of the staggered 
longitudinal susceptibility, to locate the multicritical point in the XXZ model. The 
arrow marks the asymptotic critical exponent in the Ising universality class |33j . 



two staggered susceptibilities exhibits an extremum which height is largely 
independent of system size in the vicinity of the multicritical point. Such a 
behavior seems to be consistent with a bi critical Heisenberg point. 

These observations and interpretations on the character of the multicritical 
point have to be viewed with care, having in mind controversial renormaliza- 
tion group arguments. Early renormalization group calculations in one loop 
order suggested a bicritical Heisenberg point [5|6] . Later work, in five loop or- 
der, found the biconical point to be the stable one [7f38] . however, not ruling 



out a triple point, at which at least one of the transition lines to the paramag- 
netic phase is, eventually close to the multicritical point, of first order. A most 
recent analysis indicates that all three scenarios are possible, depending on the 
model parameters [8|. In our case, A = 0.8, the Monte Carlo simulations seem 
to favor a bicritical point, with Ising and XY lines of the paramagnetic phase 
meeting the line of first order between the AF and SF phases. Of course, a 
crossover to a different multicritical behavior in the immediate vicinity of the 
multicritical point may happen. Simulational clarification of this aspect may 
need an enormous amount of computer time. 

Let us now turn to the Monte Carlo results for the full Hamiltonian, 7i{= 
'Hxxz + 7"^ca, eqs. (1) and (2), including exchange and cubic anisotropies. 
We simulated thermal properties for a few selected cases, where discretized 
biconical structures play an important role. As before, we set the exchange 
anisotropy A equal to 0.8. 

For positive cubic anisotropy, F > 0, favoring orientations of the spins along 
the diagonals of the cubic lattice, a discretized ordered BC2 phase may arise 
from the corresponding ground state, see Fig. 2. This behavior is illustrated 
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Fig. 7. Phase diagram of the XXZ model, A = 0.8, with a positive cubic anisotropy 
F at fixed field, H/J =1.8. 



in Fig. 7, when fixing the field, H/ J= 1.8, and varying F. At small values 
of F, there is an AF ordering at low temperatures. Above a critical value of 
F, F > F C (H = 1.8J) = 0.218. .J, the low-temperature phase is, indeed, of 
BC2 type. In that case, increasing the temperature at fixed F, one goes from 
the BC2 to the AF phase before then entering the disordered paramagnetic 
phase. The transitions between the two ordered phases and between the AF 
and disordered phases have been located from monitoring the staggered mag- 
netizations and susceptibilities, the specific heat, and Binder cumulants. An 
example is shown in Fig.8, where the staggered longitudinal and transverse 
magnetizations are displayed, as a function of temperature, at F/J= 1.0. 
At the transition between the AF and BC2 phases, the absolute staggered 
transverse magnetization, |m^f|, drops rather sharply, with an accompany- 
ing anomaly in the absolute staggered longitudinal magnetization, \m z st \. The 
staggered longitudinal magnetization is the order parameter of the AF phase, 
and it decreases rapidly on approach to the paramagnetic phase. 



As expected, the transition from the AF to the disordered phase seems to 
be Ising like, as follows, e.g. from the critical exponent describing the size- 
dependence of the height of the peak in xlt- The transition from the BC2 
to the AF phase is found to be consistent with the XY universality class 
[29] . Indeed, the transition between the AF and biconical phases has been 
predicted to belong to the XY universality class in an early renormalization 





0.8 




, . , | , , , , | 


iiii,i 






;I 




<|m* t |> 




0.6 




\ 

Q 
O 




_A_ 






$ 


y 


N 
X CO 

E 


0.4 




% <|m3> 


V 
V 


V 






© 


t 








$ 


Y 




0.2 




9 


y 








6 


y 








% 


y 




r\ 




*W-e- 


K- - 
~~e---e~-s& 









. . . i . . . . i 


. . . . i . 



0.5 1 

k B T/J 



1.5 



Fig. 8. Absolute staggered longitudinal and transverse magnetizations as a function 
of temperature at fixed field, H/J= 1.8, and at fixed cubic anisotropy F/J= 1.0. 
Systems with 16 3 sites have been simulated. Note that error bars are small compared 
to symbol sizes. 
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group analysis [19]. Furthermore, the cubic anisotropy is expected [33|39] to 
be an irrelevant perturbation in the three-dimensional XY case. As depicted 
in Fig. 7, there is no multicritical point of the three phases at the values of F 
we studied. 

A multicritical point may occur in the (H , T)-plane, when fixing F, F > 0. 
Indeed, our simulation results, especially at F/J = 0.15 and 1.5 [29], show a 
phase diagram comprising the AF, BC2, SF, and paramagnetic phases. How- 
ever, much more extensive simulations, being beyond the scope of the present 
study, would be needed to study the possible multicritical points in detail. Es- 
pecially, the order of the transition between the BC2 and SF phases deserves 
a careful analysis. A transition of first order, as suggested by the ground state 
analysis, would preclude a tetracritical point [19] . 



Finally, we briefly mention results for the case of negative cubic anisotropy, 
F < 0, favoring spin orientations along the cubic axes. As follows from the 
ground state considerations, see Fig. 2, biconical structures only show up at 
relatively large negative values of F. Moreover, the tilt angles usually seem to 
jump when going from the BC1 region to one of its neighboring, i.e., ferro- 
magnetic, AF, or SF, regions. 

Our Monte Carlo simulations, especially at F/J= -2 and varying the field at 
fixed temperatures [29J, suggest that there is an ordered BC1 phase which, 
however, gets destabilized at quite low temperatures. Thorough analyses of 
full phase diagrams and the types of the transition are well beyond the scope 
of the present work. 



4 Summary 



Classical Heisenberg antiferromagnets on cubic lattices with uniaxial exchange, 
A, and cubic, F, anisotropics in a magnetic field, H, have been studied, using 
ground state considerations and Monte Carlo techniques. 

In addition to antiferromagnetic (AF) and spin-flop (SF) structures non- 
collinear, biconical (BC) structures are observed, depending on the strength 
of the anisotropics, the field, and temperature. They may occur as ground 
states, which, at low temperatures, may contribute to thermal fluctuations of 
BC type or may give rise to ordered BC phases. 

Specifically, in the XXZ antiferromagnet (F= 0) biconical structures are present 
as degenerate ground states at the critical field separating the AF and SF 
states. At low temperatures near the transition between the AF and SF phases, 
BC fluctuations show up, but they do not lead to a distinct phase, in contrast, 
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for instance, to the situation for the XXZ magnet on a square lattice, where 
the AF and SF phases are separated by a narrow disordered phase. Moreover, 
in this study, the multicritical point of the three-dimensional XXZ model, at 
which the AF, SF, and paramagnetic phases meet, has been located accurately 
for A= 0.8. That point seems to have bicritical character, being the intersec- 
tion of two continuous transition lines between the AF and SF phases to the 
paramagnetic phase with a transition line of first-order between the AF and 
SF phases. Note that we did not rule out a crossover to a different character 
in the immediate vicinity of the multicritical point. 

The cubic anisotropy F may lead, depending on its sign, to a tendency of the 
spin orientations along the diagonals of the lattice, F > 0, or along the cubic 
axes, F < 0. In both cases, BC as well as SF structures have now discretized 
preferred orientations of the rry-components of the spins, thereby breaking the 
rotational symmetry of the XXZ model. 

For positive cubic anisotropy, F > 0, biconical configurations occur as ground 
states next to the SF and AF phases, giving rise to an ordered BC2 phase. 
The transition to the AF phase is found to be consistent with being in the XY 
universality class, where the cubic anisotropy is expected to be an irrelevant 
perturbation. Phase diagrams with three ordered phases, AF, SF, and BC2, 
and a paramagnetic phase have been determined. 

At small negative cubic anisotropy, F < 0, there are no biconical ground 
states. However, such structures, of BC1 type, may be stabilized at larger 
negative values of F, between the ferromagnetic and the AF or SF ground 
states. The biconical ground states then give rise to an ordered BC1 phase at 
sufficiently low temperatures. 
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